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Abstract 

We construct generalized solutions to the ultradiscrete KdV equation, 
including the so-called negative solition solutions. The method is based 
on the ultradiscretization of soliton solutions to the discrete KdV equation 
with gauge transformation. The conserved quantities of the ultradiscrete 
KdV equation are shown to be constructed in a similar way to those for 
the box-ball system. 

1 Introduction 

The ultradiscrete KdV equation: 



^ n 



(1) 



was first introduced as a dynamical equation for the so-called box-ball system 
(BBS)il][2]. Since ^ is closed on the binary values or 1, if we regard as the 
number of balls in the nth box at time step t, and assume the boundary condition 
lim = 0, it describes the time evolution of balls in a one dimensional array 

|n|— >oo 

of boxes. The updating rule of the BBS according to ([1]) can be expressed as 
follows. (See figured]) 

• Find all the pairs of a ball and an adjacent vacant box to the right of it, 
and draw arclines from the ball to the vacant box for each of these pairs. 

• Neglecting the pairs connected by arclines, find all the pairs of a ball and 
an adjacent vacant box to the right of it and draw arclines in a similar 
way. 

• Repeat the above procedure until all the balls are connected to vacant 
boxes, then move all the balls into the vacant boxes connected by arclines. 

An example of the time evolution pattern is shown in figure [51 The consecutive 
balls behave like solitons in the KdV equation. Furthermore, when we denote 
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t ...000111001001110000 



t+1 ...000000110110001110... 

Figure 1: The update rules of the BBS. A vacant box and a box with a ball are 
represented by '0' and '1', respectively. The conserved quantities are simulta- 
neously determined, {qi — 3, q2 = qs = 2.) 
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Figure 2: A solitonic time evolution pattern of the BBS. The symbol '.' denotes 
'0'. 



by qi the number of arclincs drawn in the first step of the updating rule, by q2 
one in the second step and so on, we obtain a nonincreasing integer sequence 
(qi, <Z2, qs, ■ ■ ■), which is a conserved quantity of the BBS in time[3]. 

The reason why the BBS has these solitonic properties is understood by the 
fact that H]) is obtained from the KdV equation through a limiting procedure, 
ultradiscretization[3]. The process involved in obtaining ([Ij from the KdV equa- 
tion can be explained as follows. The integrable discrete KdV equation is given 
in the bilinear form as 

(1 + S)al+\al-' = 5al-^\al+' + a^al^,, (2) 

where G 1? and (5 is a discretization parameter |_5 . By putting 

n f t-1 ' n \2 ' V'^/ 



^ gives the simultaneous equations 



{l + S)-^=Sui + ^, (4) 



(5) 
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.... -1 -2 .. -1 

-1 -2 . . -1 

-1 -2 . . -1 . . . . 

-1 -2 . . -1 . . . 

Figure 3; A time evolution pattern of a background solution. 



When we impose the boundary condition lim cr^ = const., then lim = 

|n|— )-oo |n|— >oo 

lim — 1 from dU) and ([5]) we have 

n|— >oo 

(1+5)-^ n (6) 

We introduce a parameter e and put S :— e'^^'^ . If there exists a one parameter 
family of solutions to © ""^(e), and if the limit lim elogu^(e) —: exists, 

C-J-+0 

then, using the identities for arbitrary real numbers a and b 
lim elog (e"^" + e''/') = max[a, 6], 

^hm^elog(e"/<^ •e''/^) = a + &, 

— max[— a, —b] = min[a, 6], 

we find that ?7* satisfies (H]). 

Recently, ([!]) has been studied without the constraint of binary values [51 [7]. 
One interesting feature in this extension is that there exist background solutions 
composed of so-called negative solitons travelling at speed one [5]. A typical 
example of the background solutions is shown in figure [31 Furthermore, such 
negative solitons interact nontrivially with ordinary solitons as shown in figure[4l 
A method to solve the initial value problem for this generalized situation is given 
in Ref.iTj. From these properties, the dynamical systems described by ([1]) in the 
general setting are also regarded as integrable cellular automata, just as the 
BBS. However, several important features of integrable cellular automata have 
not been clarified yet, such as the conserved quantities, the relation to integrable 
lattice models, extension to the case of a cyclic boundary condition and so on, 
all of which are well understood in the case of BBS. In this article, we report 
that a simple transformation reveals the underlying integrable structure of ([T]) 
for the generalized states. By this transformation, the relation to the BBS 
and hence to the integrable lattice models, and generalization to a different 
boundary condition become apparent. We explicitly show how to obtain the 
conserved quantities and give the expression of general solutions to ^ where 
negative solitons and ordinary solitons coexist. 
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1 2 1 ... 1 1 -1 .... -1 
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Figure 4: A time evolution pattern in which negative sohtons and ordinary 
sohtons coexist. 



2 Conserved quantities for states with negative 
solitons 

We start from the ultradiscretization of the coupled equations and By 
putting S — e^^l'^ and assuming that the limit 

lim e log <(e) [/* , lini e log <(e) 

exist, we obtain the equations 

C/*+i=min[l~C/*,F^], (7) 

= + (8) 

If we impose the boundary condition lim — lim 1/* = 0, Eqs. ([7]) and 

|n|— foo |n|— >oo 

dS]) are equivalent to ([T]). It is easily seen that Q and dH) are generalized by 
introducing a positive integer constants L and C as 

+ min [L - [/* , K!] - min [C/* , C - , (9) 
+ (10) 

Equations © and (flUl) are known as the ultradiscrete limit of the modified KdV 
equation and are interpreted as a BBS of box capacity L with carrier of capacity 
C [SI H]- At the same time, they represent the action of the combinatorial R- 
matrix of an J^^^ crystal on a certain symmetric tensor product representation as 
long as C/* G {0, 1, • • ■ , L] and Vl G {0, 1, 2, • ■ • , C}[10|. Note that Eqs. ^ and 
([To)) are closed under the values 17* G {0, 1, • • • , L) and F„* G {0, 1, 2, • • • , C}. 
Hence the time evolution patterns of the BBS are regarded as the ground states 
of a 2 dimensional integrable lattice model at temperature and its generaliza- 
tion has been investigated in detail [TTl [T^. 

When we allow U*^ or in ([7]) and (|S]) to take negative integer values and/or 
the values greater than i or C respectively, they are no longer interpreted as 
a BBS nor as a realization of the combinatorial R- matrix. For an arbitrary initial 
state including negative integer values, let M := — min V^,L — C — V^] . 

When we change variables f7* := C/* + M, V*:=V* + M, L := L + 2M, C := 
C + 2M, then we find that Eqs. ^ and ^ are unchanged under this gauge 



4 



transformation [13], that is, 



;7*+i = [/* + min [L - [7* , Vl\ - min [ul C~V^ , 

Since these equations are closed under the values G {0,1,-- - ,L + 2A/} and 
G {0, 1, 2, • ■ • , C + 2Af}, and initial values t7°, belong to the same sets, 
we can regard them as equations for a BBS with box capacity L + 2M and 
carrier capacity C + 2M, and regard them as a realization of crystal lattice 
composed of a larger tensor product representation. Thus we can conclude that 
Eqs. dH) and (ITU)) as well as Eqs. ^ and ([U do not change their mathematical 
structure even in the case where the dependent variable take negative integer 
values. 

From the above argument, we find that we can identify the ultradiscrete KdV 
equation ([T]) as a BBS with larger box capacity when its dependent variables take 
negative integer values. The conserved quantities of the BBS have been widely 
investigated in general cases and we can apply the results to this svstem fTTlfM] . 
In [TTl [TS], it is shown that the energy functions which are combinatorially 
determined, are conserved in time, and in jH] the path description for the 
conserved quantities is given. Here we give another expression of the conserved 
quantities for ([T]). 

Let M - min[[/^, V;°, 1 - [/°]. Then Eqs. © and ^ are transformed 

n 

to the BBS with box capacity 1 + 2M (with infinite carrier capacity). Let the 
number of balls in the nth box be Ul^. We consider the time evolution starting 
from the time step t — Q, and choose the two integers Mq and A'o which satisfy 
the following conditions. 

1. = Ul = for n < Mo + l imd Nq - I < n {Ul ^ fjl ^ M [n < 
A/o + 1, <'^)). 

2. iVo - Afo = (mod 2). 

3. No vacancy in the A'oth box is connected by an arcline in the first step of 
the construction of the conserved quantities which we explain below. 

The conditions 2, 3 are used to avoid the redundancy of the definition of con- 
served quantities. For the time step t, from the boundary condition, there 
exists two integers Mt, Nt such that U^ = V* = U^^ = for n < Aft — 1 and 
Nt + I < n { ul = Ul;^'^ = M {n < Mt + 1, Nt - 1 < n)). Furthermore we 
impose the following conditions on them. 

Nt'~Mt=0 (mod 2), 
Mt-Mo = t (mod 2), 
Nt-No = t (mod 2). 

At each times step t, we consider the boxes with indices from Mt to Nt- In 
other words, we neglect all the balls in the boxes n < Mt and n > Nt. Then, 
as in the case of the BBS with box capacity one, we draw arclines from balls to 
vacancies to determine the values of the conserved quantities. (See figure O) 
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• Draw an arcline from one of the leftmost balls to one of its nearest right 
vacancies. 

• Draw an arcline from one of the leftmost balls which are located on the 
right of the vacancy we have drawn the arcline to one of its nearest right 
vacancies. 

• Repeat the above procedure until one has drawn arclines between all the 
pairs of balls and vacancies. Note that if we drew an arcline to a vacancy 
at nth box, then the ball from which we draw an arcline in the next step 
is located in the box whose index is greater than n, and that there must 
be no arcline starting from a ball in the last box due to the condition 3. 
Let qi(t) be the number of arclines drawn in this step. 

• Then neglecting all the boxes and vacancies which are connected by ar- 
clines, repeat the above procedure and let q2{t) be the number of arclines 
drawn in the second step. 

• Repeat the above procedure until all the balls but those in the A'tth box 
are connected and let qj (t) [j = 3,4,...) be the number of arclines drawn 
in the jth steps. 

• Let qj := q,{t) ~ ^{Nt - Mt) [j = l,2,...,2Af) and q, qj{t) {j = 
2M+ l,2M + 2,...). 

Note that we obtain the updated state of the BBS with carrier capacity C if we 

exchange the balls with the vacancies connected to them up to the Cth steps, 

c 

i.e., if we exchange the first 9 j pairs of balls and vacancies. Then we find 
the following proposition. 
Proposition 1 

The sequence fgi, q2, Qa, . • ■) does not depend on t. 

The statement of the proposition is the direct consequence of the fact that 

k 

qi is equal to the number of balls moved by the carrier with capacity k at 

1=1 

one time step which is the value of fcth energy function [11] [12]. Here we map 
the system to the BBS with box capacity 2M + 1 in which the number of balls 
is finite, that is, the nth box is vacant for n < Mt — 1 and Nt < n. Of course 
we can make another choice. We could impose the condition Nt — Mt = 1 
(mod 2) or the condition that a vacancy in the iVoth box is connected by an 
arcline in the first step. This redundancy is caused by the fact that there are 
essentially four types of finite systems which are equivalent to this system. Once 
we determine the type however, the sequence ((71,52, • • ■) does not depend on 
Nt and Mt as long as they are sufficiently large and, when there is no negative 
soliton q2M+j — qj {j = 1,2,.. .), where qj is the jth conserved quantity defined 
in the original BBS. 
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Figure 5: The scheme to count the conserved quantities of the BBS with larger 
capacities. We have to choose different sequences of boxes at time step t and 
t + 1. From this figure, we have qi{t) = 7, 52 (^) = 5, 53 (t) — 2, q4{t) = 
1, g5(i) = l, <76(i) = landgi(t+l) = 7, 92(^+1) = 5, 93(^+1) = 2, 94(^ + 1) = 
1, q5it+l) = f, qeit+l) = 1. Since M = 1, Nt-Mt = f4 and A^t+i-A//t+i = 14, 
we obtain qi — 7 — ^ = 0, 92 = 5 — -y = —2, (73 =2, (74 = 95 = = 1 in both 
time steps. 



3 Soliton solutions interacting with a background 
state 

For the original BBS where every box capacity is one and only a finite number of 
balls exist, all the solutions are proved to be constructed by ultradiscretization 
from the soliton solutions of the discrete KdV equation ©[HI- Although the 
ultradiscrete KdV equation ^ over arbitrary integers is shown to be equivalent 
to the BBS with larger box capacity, the solutions to them can not be the 
ultradiscrete limit of soliton solutions to (O, because the boundary condition 
of the BBS changes such that the number of balls becomes M for |n| 00. 
Nevertheless we shall show that solutions can be constructed from the soliton 
solutions to © through ultradiscretization with an appropriate scaling limit. 
Here we treate only the case M = 1. Our method for constructing the solutions 
is equally applicable to the case M > 2, although the analysis of that case is 



more elaborate. 

The coupled equations for the BBS are explicitly written as 

C/*+i-min[3-C/*,V;*], (11) 
V^+, = Ul-, + V,\-U'+\ (12) 

with the boundary condition 

hm C/*- lim V^, = l. (13) 

71 1 — ^ oo I n I — ^ oo 

From and (O, if we put 5 — e^^l"^ [L eX^) and assume that 

£-S-+0 

exists, we find 

P'n+i + P'n^ = max + -L,p\ + . (14) 

From ([3]), if we put 



C^n ^Pn-Pn - Pn+1 + (15) 



they satisfy PT|) and (jl2p for L = 3. Hence we will construct solutions of 
(fT4)l which satisfy the boundary conditions (IT3l) . 

First let us consider a state where only negative solitons exist. The corre- 
sponding solution is called a background solution [7] . It is a stationary solution 
of the form C/^ = F{n — t) where {F{ky\^^_^ is a 01 sequence with only finite 
number of Os. Let the number of Os be s + 1 (s G Z>o), and let 60 = 1^0 = 0, 
h e Z>o {i = 1,2, ...,s) and 5i G {0,1} [i = 1,27- ..,s). We define a 01 
sequence {F[k)}^^_^ as 

^(^) ^ ( k = - ELo(2&. + '5. + «) (j = 0, 1, 2, ... , s), 
\ 1 otherwise. 

Apparently any background solution F{k) can be expressed as F{k) — F(k — a) 
for some integers a, s, {6^} and {Si}. 

The iV soliton solution to ([2]) is given as 

(18) 

where pi G C and 7^ G C (i = 1, 2, . . . , n) are arbitrary but pi ^ pj for i ^ 
j[4l [5]. We often denote [fc] := {1, 2, 3, . . . , fc} for a positive integer k. To 
take the ultradiscrete limit, let b = e^^^*^, pi ~ NiC^^'^'^^ ~% — MiC^^^^"^ 
(i = 1,2,...,N) where C, G K, L, P^, A^^, Af,; > and A, ^ Nj, M, ^ Af, for 
i ^ j. Then taking the limit lim elogcr^ —: pl^, we have 



p„ = max < 0, max 



^(a + t/^, - nmin[P„ L]) - ^ min[P„ P,] 



(19) 
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For i = 1, we find that C/* and V* satisfy © and ©, and hence ([T]). In 
particular, if Pi G Z+, G Z, ([T^ gives an N soliton solution of the original 
BBS. 

We are, however, concerned with Eqs. (|TT]) and (|T^ with the boundary con- 
dition (jl3p . Hereafter we put L — 3 and prepare several symbols and notations: 

/:=^(5,, = s - ^ + 1 (mod 2), I' -.^ I + Ss+i, 



K 



j=0 i=o 



We also assume that Si — I for i — mi, m2, . . . ,1111 (1 < mi < TO2 < • • • < to; < 
s), and m;+i := s + 1. Let us consider an rioo soliton solution (rioo ;3> 1) of the 
form ([H]). We assume that Pi = 2 {1 < i < - I') and i^' := Pi+n^-i' = 1 
(1 < i < I'). The phases are chosen as 

r 1 (1 < i < no), 

-2j + l (no + S,-i + l <f <no + B, (j = l,2,...,s)), 
[ -2s- 1 (no + B, + l<i), 

C'i = d+n^-i' ■= -{mi -i) {i = 1, 2, /')■ 
Then we find 



max 

0<fc2 



fcl 



max 

0<fc2 
fe2 



Y^{C.-2{n-t))+Y,{C[-{n-t)) 

4=1 1=1 

2fci(fci-l)-/s2(fc2-l)-2fcifc2 

-(2fci + fc2)(n - t + 2no) +51'^' 



q - 2fci(fci - 1) - fc2(fc2 - 1) - 2fcifc2 

i=l 

— 2no(n — t + 2no) + 2riQ + 3no) 



where 




(fc < 0). 



It is obvious that if is a solution to (jl4p . ^^^"5, + at + /(n) is also a solution 
for an arbitrary function f{n) and arbitrary values a, b and c. Hence we find 
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that p*j — Prio.Moo 

{n - i) is a solution to dH]), where 



0<fc2 

fel ' fc2 

-(2fci + fc2)a: + ^C, + 



i=l 1=1 



2ki{ki - 1) - /C2(A:2 - 1) - 2/ci/c2 



Since no, noo are large but arbitrary positive integers, the limit 



Pf{x) ■■= lim lim Pno,n^{x) 



max 

ki GZ 
0<k2<l' 



-(2/ci + A:2)x + ^Q 

^ q - 2fci(fci - 1) - fc2(fc2 - 1) - 2fcifc2 



fc2 



(20) 



is also a solution to (fTH) . 

When we rewrite in ([T?l) as 



-ra; + max 

0<k2 <l' 



4=1 



and using the relations such as 

flo + 02 > 2ai — > max[— 2a; + a2, — a; + fli, oo] = max[— 2a; + 02, flo] 

Eq. (PP)) can be simplified as shown in the Proposition [5] below. First let us 
define some symbols and notations. 

s I 

Ct ■■= Y.(2j - 1)6, + ^(m, - j), Qo := Ct + if (if - 1) + (if - 1){K - I - I) 
5i {mi,TO2, ■ • ■ 5o:=[/]\S'i, I:=[2K-l]\I 



/ := |m e pif- 



m = 2B,_i + A,_i + 2k, + 1, fc» = 0, 1, . . . , 6,_i - 1, i G 5o 



The quantity 8^ is defined for r G i as 

• If i e 5o and 2Bi_i + Ai_i + 2 < r < 25.; + A^, then putting r 
2B,_i + A,_i +2fc (fcGZ), 



e.;=5] (2j-l)6, + (j--^<5^)5, 



+ {2i - l)k 



(B,_i +k + A,_i)(B,_i + fc + A,_i - 1) + (B,_i + fc)(B,_i + fc - 1). 

(21) 
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• Hie Si, 2Bj_i + Ai_i + 1 < r < + A, and r = 2B,_i + A,_i + 2k 
[k e Z), then is given by (1^ . 

• If i e 51, 2B,_i+A,_i + l <r < 2Bi + A, and r = 2Si_i + Ai_i + 2fc + 1 
(fc € Z), then 

i-l i j 

Or ■■= - 1)6, + - E + (2* - 

+ + fc + A,)(S,_i + fc + A, - 1) + + fc)(B,_i + fc - 1). 

(22) 



Proposition 2 

PF(a;) := max p^f\x) , p''p\x) , p^p\x) 

where 

p^f\x) :=max[(2x- l)fc - 2fc(fc + 1)] , 



(23) 



(24) 



max [-{2K - I + 2q)x - q{2q + AK - 21 + 2s - 1) - Qo] {Ss+i = 

q>l 

max \-{2K - I + q)x - {q - 1) + 2K - I + s) - s - Qq] {S,+ 
9>i L \2 / J 



max 

l<r<2K-l 



p^^\x) 



(25) 
(26) 



and the function p^"^^ (x) is given as 



rx — Qr r e I, 
ril. 



From (fTS)) . if we define 

Uf{x) := pf{x + 2) - 2pf{x + 1) + pf(.t), (27) 

17* = UF(n — t) gives a solution to pT|) and ([H]). By using the fact that 
Pf{x) is a continuous and convex piecewise hnear function, and by calculating 
the intersection points of each linear functions appeared as linear pieces in the 
function Pf{x), we obtain the explicit form of Uf{x). Let us define ti{x) {i = 
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1,2,...) as 



t2(x) 





{-2<x< -1), 




(-1 < X < 0), 


\ 


otherwise. 


[ 2x + 


5 (-5/2 < a; < 


I -2x- 


-1 (-3/2 < a; < 


l 


otherwise. 


2k 








1=0 




k-1 




Y^hix- 


f 2i). 


i=0 





By straightforward but a Uttle tedious calculation, we find 

Uf{x) = U^p\x) + f/f (x) + 4+^(x), 



where 



k=l 



J^Mx + K + s + 2k) ((5,,+i=0), 

oo 

Y,tiix + K + s + k) ((5,+i = l), 



fc=0 



Wp'ix) := *2&.+5. + 2Sz-i + A,_i +1-1). 



(28) 



(29) 



(30) 



(31) 



An example of the function Uf{x) is shown in figure [51 Comparing these ex- 
pressions with F{k), we find that UF{k) — F{k) for k E Z. Thus we have proved 
the following proposition. 

Proposition 3 

It holds that UF{k) = F{k) for k e Z. 



Furthermore, comparing F{k) with the conserved quantities {qj} defined in 
the previous section, we have the following proposition. 

Proposition 4 

The number of negative soUtons is equal to —qi — q2, i.e., s + 1 = —qi — q2, and 
qi — 92 is the number of solitons with amplitude 1 which is used to construct 
the background solution, i.e., qi ~ q2 — I' ■ 



Next we consider the solutions in which m ordinary solitons exist in a back- 
ground state. Since N soliton solutions are given by (ITOl) and the parameters 
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Up(x) 




F(x)= 1111011110111011111111 

Figure 6: An example of F{x) and Uf{x). The parameters are 61 = 1, 5i 
1, b2 = 2, 52 = 0. 



for the background solution are already known, we have only to add new pa- 
rameters {Pi}^i C^i Pi > 3) and {Ci}^^. When we denote the set of indices 
for a background solution B, then the solution is given by 



p„ = max 

JiC[m],J2CB 



J2 (- min[3, P^]n + Prt + C,) + ^ i-P^in - t) + C,) 



.ieJi 



min[A,P,]- J2 niin[P„P,]-2 ^ min[Pj, Pj] 

i¥^j:i,j&Jl i^j,i-j&J2 i&Jl,j&J2 



max 

JiC[m],J2CB 



Y (-3n + P,t + C,) + Y, (-^^(" + C^) 



.ieJi 



ie-h 



J2 min[P„ P,] - J2 ^MP^,Pl] - 2 51 1-^1 \Pi 



j'e Ji 



je.h 



max 

JiC[m] 



^ (-3n + P,t + a) - J2 ' 



+ max 

J2Ci3 



ie./2 j#i,«,ieJ2 



Here we used the abbreviation 



max [• • • 1 = max < 0, max [■ • • 1 

J<ZA 1 J<ZA,J=^!ti 



Therefore we obtain the following theorem. 



Theorem 1 

The solution to (jl4p representing the interaction of m soliton solutions with the 
background given by ((23)) . has the expression: 



p„ — max 

JC[ml 



5]](-3n + P^t + C.,) - Y min[A, Pj] + - t + 2\J\ 



ie.j 



(32) 
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..22222.0. .22.0. ...0 

232. . .220. ...0 

0.222. .32. ...0 

0... 220. 222220 

0....30 322 

0...03 0.. 22222.. 

Figure 7: A solution [/* given by Here {/* = - p\-^ - p^+i + p*,+\. The 
symbol '.' denotes '1'. The parameters are Pi =4, Ci = —7, P2 — 1, C2 ^ 
—5, bi = 2, Si — 0, 62 = 1, S2 — 1. When we make gage transformation as 
?7* — 1 — > l/^, this time evolution pattern coincides with that of figure SI 



The multi-soliton solutions with a background is also given by Y. Nakata in an 
algebraic form (nested maximization) with the ultaradiscrete vertex operators [TB] 
So far we do not find the exact correspondence between our solution p2p and 
his formula. 

From (1321) , we find that the phase shift of negative solitons after collision of 
an ordinary soliton is —2. An example of a time evolution pattern is shown in 
figure [71 Since the solution ((5^ shows that the ordinary solitons move freely 
for t 3> 1 as those in the BBS with box capacity one, comparing with the 
construction of conserved quantities, we find the following proDOsition[l7). 

Proposition 5 

For a solution to (H]) determined by (15^ . let Y be the Young diagram corre- 
sponding to the partition (ga, 94, 95, • • • ) constituted by the conserved quantities, 
that is, the length of its ith column ofY is qi+2- Then, the length of jth row 
in Y is equal to Pj — 2. 



4 Concluding remarks 

We constructed the conserved quantities of the ultradiscrete KdV equation ([T]) 
with negative solitons by using the transformation to the BBS with larger box 
capacities. The solutions to (H]) in which negative solitons and ordinary solitons 
coexist are also obtained as the solutions to this extended BBS. In the BBS 
with box capacity one, every state is obtained by ultradiscretization of a soliton 
solution of the discrete KdV equation, and using this fact, initial value problem 
is solved with elementary combinatorial methods. We conjecture that this fact 
also holds for the extended BBS and can be used to solve the initial value 
problems. Investigating this conjecture is a problem we wish to address in the 
future. 
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